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ABSTRACT 

The interaction of a 1 kW arcjet thruster plume with a communications sig- 
nal is evaluated. A two parameter, source flow equation has been used to 
represent the far flow field distributions the arcjet plume in a realistic 
spacecraft configuration. Modelling the plume as a plasma slab, the interac- 
tion of the plume with a 4 GHz communications signal is then evaluated in terms 
of signal attenuation and phase shift between transmitting and receiving 
antennae. Except for propagation paths which pass very near the arcjet source, 
the impacts to transmission appear to be negligible. The dominant signal loss 
mechanism is refraction of the beam rather than absorption losses due to colli- 
sions. However, significant reflection of the signal at the sharp vacuum- 
plasma boundary may also occur for propagation paths which pass near the 
source. 


INTRODUCTION 

A low power arcjet thruster is currently being developed as a high per- 
formance alternative for auxiliary propulsion applications such as north-south 
stationkeeping on communications satellites (refs. 1 and 2). When operated 
with storable propellants at power levels between 1.0 and 1.5 kW, the arcjet 
can provide specific impulse values greater than 400 sec, offering potentially 
large benefits in terms of spacecraft mission life or payload capability. A 
laboratory arcjet thruster has demonstrated technology readiness in terms of 
lifetime and performance (ref. 1). 

As the arcjet moves closer to being a flight-ready propulsion system, 
spacecraft integration concerns become of significant interest. Many of the 
key issues center around the arcjet plume, because of the potential impacts 
that the plume may have on other spacecraft subsystems. Potential impacts of 
primary concern include: electromagnetic contamination due to radiated emis- 

sions from the arc discharge and plume, induced environment effects on communi- 
cations signals, thermal loading of sensitive spacecraft surfaces, and momentum 
losses due to plume, impingement. The relative importance of the above concerns 
depends upon the actual thruster-spacecraft architecture. 

In order to accurately assess the arcjet plume impacts, a flow field char- 
acterization of the exhaust is required. A partial description of the charged 
fraction of the arcjet plume has been obtained using electrostatic probes 
(ref. 3). Measurements 30 cm downstream of the thruster exit plane indicated 
plasma number densities on, the order of 5xl0^/cm3 at electron temperatures less 
than 1.0 eV. The exhaust is only si ightly. ionized (<1 percent). An increase 
in measured plasma .number density was observed with increasing arc current at 
a constant mass flow rate and also with an increasing mass flow rate at con- 
stant current. Radial measurements 18 cm downstream of the thruster exit plane 



indicated that the plasma number density decreased exponential ly off centerline 
while values of electron temperature remained relatively constant with radial 
position. Similarly, an axial survey of the plume centerline revealed that the 
plasma density decreased as the inverse - of the distance squared. 

This paper focuses on potential communications impacts due to the induced 
environment of an arcjet plume. Due to physical limitations of a spacecraft 
structure, it may be necessary for uplink and/or downlink common carrier sig- 
nals to be transmitted through a portion of the plume, as illustrated in 
figure 1. The interaction of the radio frequency (RF) signal and the charged 
species of the plume must be well understood.. Typical RF signal-plasma inter- 
actions include: reflection of the transmitted signal at the vacuum-plasma 

boundary, attenuation and phase shift, and generated noise on both signal 
amplitude and phase. 

Potential arcjet plume impacts are difficult to quantify. However, the 
interaction between a plasma exhaust and common carrier signals has been stud- 
ied previously (refs. 4 and 5). During the Solar Electric Propulsion develop- 
ment program, the transmission of 2.2 GHz signals through a 30 cm mercury ion 
thruster exhaustwas investigated both analytically and experimentally. It was 
reported that signal amplitude loss and phase advance increased with ion beam 
current.' Maximum transmission loss was measured at 0.38 dB or 8 percent. 

Phase shifts advanced as high as 67°, although the effect of this shift was 
only significant during fast ion beam transients. No additional noise was 
observed on the carrier frequency during normal thruster operation, and the 
carrier frequency was not offset due to plasma effects. Poor space simulation 
may have existed in the experiments due to exposed metal vacuum tank surfaces 
causing signal reflections and an abnormal EM environment. Further, it was 
concluded that improvements in the analytical modelling of, the plume electron 
density distributions and RF beam spreading were needed to confidently predict 
RF signal-plume interactions. 

This paper presents a first order evaluation of the potential impacts of 
the arcjet thruster plume on communications. Assuming that gas dynami c expan- 
sion, rather than recombination, predominates in the far flow field plume, a 
two parameter equation (ref. 6) may be used to mathemati cal ly represent the 
electron number density distribution in the arcjet exhaust. Available experi- 
mental data are used to determine the values of the plume parameters.. The 
interaction of the plume with a communications signal may then be evaluated by 
treating the plume as a homogeneous, isotropic plasma which perturbs the exist- 
ing electromagnetic field state in space. .The two primary properties which 
affect the steady state interaction are electron number density and tempera- 
ture. Other factors include the RF power density and frequency, the propaga- 
tion path (spacecraft geometry), static magnetic field intensity, and ion 
density and mass . 

The first portion of this paper summarizes the basic theory for the inter- 
action between a plasma and a radio, frequency signal . The rationale for and 
the application of a source flow, two parameter equation to mathematically 
represent the arcjet plume is al so presented. Comparisons are made between the 
analytical’ model and experimentalprofi les. Using a realistic space flight 
geometry, the communications impact is then evaluated for a 4 GHz signal trans- 
mitted along propagation paths running parallel to and intersecting the plume 
centerline. Transmission attenuation and phase shift are, presented as a func- 
tion of propagation path. 
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NOMENCLATURE. 1 


A location of antenna .behind ar.cjet nozzle exit plane, 1 m 

A(N,\) plume normalization constant (eq-., (10)): 

B distance of antenna away from plume axis.% m 

c speed of light, m/sec 

D assumed plasma slab dimension, m 

D a antenna aperture dimension, m 

Dp Langmuir probe diameter, m 

e fundamental charge, C 

f (9) functional distribution of density ■> 

I current, A 

L distance from antenna aperture to vertex, m ■ 

m charged specie particle mass, kg 

M simplifying relation for (1/p - 1) 

n number density, particles/m3 

h source number density, particles/m^ 

N plume parameter 

n c critical electron number density, electrons/m^ 

P 2 x the entrance ordinate (geometrical optics approx.), 

r radial distance from arcjet source, m 

R radius from plume centerline, m 

distance from antenna aperture to center of plasma, m 
R e exit radius of nozzle, m 

T transmission loss, dB - . 

v velocity, m/sec -v.«- • - • • • 

V voltage, V : : 1 .••••: *vv 

Vp probe voltage, V ^ ■*.; 



X 


distance along transmission path'-, . m 
x 0 intersection of transmission path with plume boundary, m 
z axial distance along plume centerline, m 

9 angle off plume centerl ine,- degrees 

e 0 permittivity of free space, F/m 

X plume parameter c 

X s signal wavelength, m 

vi plasma index of refraction 

v collision frequency, rads/sec 

<P o angle of propagation path with respect to centerline, deg 
A4> phase shift of RF signal, deg 

a) signal frequency, rads/sec 

cop plasma frequency, rads/sec 

up gyrofrequency, rads/sec , . • 

Subscripts or Superscripts: . *.• : " 


a 

arc 

e 

electron 

i 

ion 

1 im 

1 imiting 

max 

maximum 

0 

centerl ine 

P 

plasma 


ANALYSES 

RF Signal-Plasma Interaction . .■■ ■■■■ 

Free electrons in the plume, because they are lightweight and very mobile, 
are the major source of interaction with RF transmissions as they readily 
absorb more energy than ions. However, at low frequencies, particularly below 
500 MHz, the ion contribution also becomes important. RF transmission loss and 
phase shift result from a combination of physical processes including reflec- 
tion of the signal at the vacuum-plume boundary, absorption, and refraction. 
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Reflection of the transmitted signal by the plume depends upon the dielectric 
characteristics of the plume boundary. Generally, if the plume characteristics 
vary gradually near the boundary, the plume may be assumed to be adiabatic and 
reflection losses are minimal. Absorption, or 1 ine-of-sight attenuation, 
removes energy from the signal while refraction (beam spreading) merely 
redirects the incident energy. 

A detailed theoretical treatment of the interaction between a plasma and 
microwave transmission signals is given by Heald and Wharton (ref. 7). The 
following analysis simplifies the interaction through the use of several criti- 
cal plasma assumptions. A complete discussion of the limitations and validity 
of these assumptions will be presented later. As described by reference 7, the 
plume may be treated as a homogeneous, isotropic plasma which perturbs the 
existing electromagnetic field state in space. In a uniform medium without a 
magnetic field, a complex dielectric constant may be used to describe the prop- 
agation of plane electromagnetic waves. For a signal frequency u, a critical 
number density n c is defined by 



e mu 
0 



(1) 


For n < n c , the plasma behaves as a nearly transparent dielectric; for 
n > n c , the medium is totally reflecting. When the signal frequency is greater 
than the plasma frequency, the plasma is a relatively low-loss dielectric. In- 
the limit where the square of the effective col 1 i sion frequency v 2 << u 2 -u p 2 
and v 2 << u 2 (u 2 -up 2 ) 2 /up 4 , the plasma index of refraction (real) is 



( 2 ) 


where the index of refraction for a perfect vacuum is 1.0. 


Using a simple adiabatic analysis of a plasma slab, the phase shift caused 
by the presence of the plasma along a transmission path is 


A4> 


r 00 

r n (x)l 

1- 

1 - e 


n 

L c J 



(3) 


where the propagation path extends from the transmitting to the receiving 
antenna. The shift is a phase advance since the index of refraction of the 
plasma is smaller than unity. For n << n c , the shift is linearly proportional 
to the electron number density averaged along the propagation path. 

Assuming that the plasma is at a distance of at least one wavelength away 
from the antenna to avoid induction effects and the dimensions of the plasma 
are large in comparison to a wavelength, the geometrical optics limit is an 
appropriate approximation fora two-dimensional analysis. The geometry is 
illustrated in figure 2. The plasma is assumed to be a uniform cylinder of 
diameter D which constitutes a divergent cylindrical lens. The distance P/2 
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is the entrance ordinate for the highest angle ray refracted by the plasma 
which would be received by. the antenna. The- transmission: loss -"due to refrac- 
tion (assuming negligible absorption and reflection) , i.s written in terms, of the 
ratio of power received by the antenna with and without the presence of a 
plasma. Specifically, the transmission ratio may, be written as 

T - £ °<L_L,za) (4) 

P<H = 1) 4M«(L , St > + D(L + 2 St) 

where M = 1/p - 1. For P/D << 1 , the cylinder is equivalent to a slab of 
thickness D. It is conventional to describe transmission loss in units of 
decibels where T in dB = 10 log (transmission power ratio). 

As mentioned previously,, for u >> up, absorption losses in the plasma may 
be assumed negligible. However, for comparison purposes, this transmission 
loss may also be calculated, in dB, using 

T' dB - -4-3 ^(j 5 ) (5) 


Similarly, an estimate of the transmission loss due to reflection at the 
sharp vacuum-plasma boundary may be obtained by 







where P/D << 1 is the equivalent of a plasma slab. 


( 6 ) 


Flow Field Description of the Arcjet Plume 

A detailed numerical treatment of plume expansion requires an accurate 
description of the nozzle flow field as an input. The arcjet thruster is char- 
acterized by a nonequilibrium, highly viscous, compressible flow field where a 
considerable amount of heat addition occurs in the nozzle. Although computa- 
tional methods exist to analyze the nozzle flow of higher thrust (> 22 N) pro- 
pulsion devices, the complex nature of arcjet gas/plasma dynamics and the small 
nozzle size precludes the immediate adaptation of conventional numerical tech- 
niques. In order to obtain reasonable predictions of the arcjet plume impacts, 
it is helpful to begin with source flow analyses. 

In the far flow field, gas dynamic expansion, rather than recombination 
processes, dominates the plume expansion (ref. 8). The plasma density in the 
plume may then be treated as though it originated from a point source. Thus, 

"p v (V) ‘ II n p <r;e> '' dA ' <7) 

where np is the source, plasma density, v is the species velocity, and 
R e is the exit radius of the. arcjet- nozzle ; . .Because recombination of species 
in the near. field of the arcjet source is significant, the "source". density 
defined here is considered a virtual , source . 

It is assumed that the density at any location in the plume n p (r,6) can 
be related to its centerline value n p (r,0 = 0°) by a function f(9). The two 



parameter, functional distribution of Reynolds (ref. 6), converted to spherical 
coordinates, is used, yielding the relation ■ 


n D <r,e> 

n p (r,6 = 0°) 


= f (0) = e 


-{[X(l - cos 0) ] n } 


COS 0 


( 8 ) 


The values of the parameters N and X are determined experimentally. 
Substitution of equation (8) into (7) gives the integral form. 


n p V 


(* R e) - n p (r,0 = 0°) 


2ir 0 


u 


^ vf (0)r^ sin 0 d0 dtp 


(9) 


oo 


where 0] ^ m is the limiting angle to the flow of charged species, 


Since the free electrons in the plume are the primary interaction mechan- 
ism, the subsequent analysis will consider only the electron number density 
distribution. As the electrons lose very little energy during expansion, the 
electron velocity remains relatively constant. Defining a plume normalization 
constant A(N,X) as 


A(N,X) = 1 



sin© f (0)d0 


(10) 


the plume structure, based on this far field approximation, is described by 

2 


n e <r,0) 


= A(N,X) fpj f(0) 


( 11 ) 


Figure 
function of 
and 1 . For 
describes a 


3 presents normalized profiles of electron number density as a 
the angle off centerline for several values of X where N = 0.8 
a given value of either N or X, increasing the other parameter 
more confined plume profile. 


Estimation Of Arcjet Plume Communications Impacts 

Communications satellites are generally located in geostationary orbi ts so 
that their position relative to a viewer on earth remains constant. Typical 
communication frequency pairs are at 4 and 6 GHz and at 12 and 14 GHz. The 
lower value of each pair is always used as the downlink transmission frequency 
because the absorption characteristics of the atmosphere increase with fre- 
quency. It is noteworthy that, while the atmosphere absorption impacts 
increase with frequency, the effect is opposite with a plasma. This is due to 
the fact that the refractive index for a plasma,' in contrast to typical dielec- 
trics, is less than unity. As shown in equations (1) and (2). for a given n e , 
the index of refraction approaches unity as n c (or rather, uO increases. A 
4 GHz signal (X s = 7.5 cm ) was used in this analysis as it represents the 
worst case . ’ 


Typically, core structures for communications satellites have dimensions 
on the order of a few meters. Figure 4 illustrates the thruster/antenna con 
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figuration used in this analysis. The distance A corresponds to the location 
of the antenna behind the thruster exit plane while B denotes the distance away 
from plume centerline. The angle <p 0 denotes the angle of the propagation 
path with respect to centerline, where cp Q = 0° corresponds to a line of sight 
parallel to the plume centerline and cp Q = +76° corresponds to intersection of 

the RF signal and the plasma source. The antenna and the thruster lie in the 
same plane. Calculations were made for propagation paths running both paral- 
lel to and intersecting the plume centerline. For the first case, A and B 
were constant and the angle <p Q was varied; in the second case, A was con- 
stant and B was allowed to vary away from thruster centerline with <p = 0° . 

The parameters listed in table I were used to define the representative 
arcjet plume for the two cases. Using equation (11), the electron number den- 
sity was then calculated at discrete intervals along a given propagation path. 
The phase shift was obtained by direct numerical integration of equation (3). 
Although the electron number density varied significantly along a given propa- 
gation path, the transmission attenuation was calculated by treating the plume 
as a plasma slab of uniform density. The relevant plasma dimension (diameter 
or slab thickness) was determined by calculating the distance over which the 
density decreased by a factor of 10. The sensitivity of the calculations to 
the selection of plasma dimension will be presented later. 

It was desired to investigate only the RF signal-plasma interaction and 
not antenna-specific effects that become important during actual thruster inte- 
gration. Therefore, the antenna size and placement selected for the analyses 
were such, that it could be assumed that the interaction occurred in the antenna 
far field, that there were no induction effects, and that the entire RF beam 
intersected the plasma. Also, it was assumed that the RF power was low enough 
so that a linear variation of power density existed in the plasma; i.e., the RF 
signal did not perturb the plasma medium. A list of the major assumptions used 
in the plasma model and the calculations to follow may be found in table II, 
along with the criteria necessary to prove the validity of each assumption. 


APPARATUS AND PROCEDURE 

The Langmuir probe is a well established technique for the measurement of 
plasma properties (refs. 9 to 11). A previous investigation (ref. 3) has 
explored the application of the Langmuir probe to the measurement of plasma 
number density and temperature in the arcjet plume. Briefly, the Langmuir 
probe is simply a conductive element of known geometry which is -connected to a 
power supply capable of biasing the probe to voltages both positive and nega- 
tive with respect to the plasma to be investigated. The current drawn to the 
probe as a function of applied voltage yields fundamental information on the 
plasma properties, primarily electron number density and temperature. 

A spherical Langmuir probe with a diameter of 6.35 mm was used to obtain 
radial profiles of electron number density at various distances downstream of 
the arcjet nozzle exit plane. It has been shown (ref. 3) that the ratio of 
electron number density to the centerline value may be obtained with reason- 
able accuracy by monitering the electron current collected by probe at a 
voltage positive with respect to plasma potential. To avoid the necessity of 
taking discrete V-I characteristics at each measurement position, this method 
was used to obtain the profiles. A complete V-I characteristic was obtained 
on centerline at z = 18.4 cm. 
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The arcjet thruster was a laboratory model (ref. 3), operating on a 1:2 
mixture of nitrogen to hydrogen to simulate fully decomposed Hydrazine. For- 
the data presented, the arcjet thruster operating conditions are given in 
table III. The arcjet anode was connected to facility ground. 


RESULTS- AND DISCUSSION 

This section presents a comparison between measured electron number 
density profiles and the source flow analysis. An explanation for the unusual 
distribution function of plasma number density is proposed. The interaction 
of a 4 GHz signal transmitted through a plasma with properties characteristic 
of an arcjet plume is then analyzed by modelling the arcjet plume as a plasma 
slab. The relative importance of the dominant loss mechanisms are compared, 
and illustrations of the number density profiles that selected propagation 
paths would encounter are presented. Finally, an evaluation of the transmis- 
sion impacts is provided in terms of signal attenuation and phase shift. The 
impacts are presented as a function of propagation path for signals running 
both parallel to and intersecting the plume centerline. 


Comparison of Experimental Profiles with Two Parameter Equation 

Figure 5 compares calculated number density profiles, normalized to the 
centerline value, with experimental data. The profile data were obtained by 
measuring the electron current collected by the Langmuir probe at an applied 
voltage of +3 V. Measurements were made at various distances downstream of 
the thruster exit plane. Table III lists the arcjet operating conditions and 
plasma plume parameters for the comparison. A complete V-I characteristic 
obtained at the centerline position for. z = 18.4 cm yielded an electron number 
density of 7xl0^/cm3. The values of the plume parameters which best describe 
the profile data taken at z = 18.4 and 28.6 cm are N =0.8 and X = 4. With 
appropriate selection of N and X, there is very good correspondence between 
the experimental, data and the analytical approximation in the far field. In 
the near field, the centerline values also agree reasonably well, as shown in 
figure 5(c), but the parametric equation greatly overpredicts the electron 
number density along a radius when the axial distance z is less than 10 cm. 
The values of N and X which would provide a best fit of the data at 

z = 8.3 cm are 0.9 and 7, respectively. 

An estimate of the background plasma number density is also shown in 
figure 5 as lO^/cm^. The measured number density reached the background plasma 
density at radial distances corresponding to an angle of about 75°. A complete 
vacuum expansion would probably show a decrease in density beyond 75°. Meas- 
urements of the arcjet plume expansion in a facility of higher vacuum will pro- 
vide a more accurate description of the high angle variation in plasma density. 

Unlike traditional source flow descriptions of plumes, the arcjet plume is 

characterized by a slowly decaying,, exponential decrease in electron number 
density off centerline. In addition, the heavy particle -ions follow -the same 
radial pattern as the electrons. v This is in direct contrast to the highly con- 
fined plumes of chemical rockets or ion thrusters where a cosine distribution 
would provide a better, descri ption of the angular variation of density (refs. 6 
and 12). A cosine distributiqn is equivalent to values of X >20 in the func- 
tion f(9) given by equation (8). This latter situation is illustrated in 
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figure 3(a) using X = 100.. However, cumulative arcjet plume data have indi- 
cated values of plume parameters ranging from N ;== 0.7 to 1 .0 and X = 3 to 7. 
There are two possible reasons for the observed gradual decrease in plasma 
number density off centerline: (1) facility effects are dominating the expan- 

sion process; or (2) the functional description is real and valid and is a 
result of the complex, highly viscous flow generated in the arcjet nozzle. 

Both possibilities deserve comment. 

Firstly, although the facility pressure is relatively high ( 2 xl 0 _/ * Torr), 
the location of the first mach disk is greater than 100 cm downstream of the 
thruster exit plane. It does not seem likely that the background pressure is 
affecting measurements within 65° of centerline and certainly not within 40° 
where a cosine distribution would be readily apparent. Secondly, it seems 
unlikely that inelastic collision processes would alter the charged specie dis- 
tribution since the collision frequencies for excitation and charge exchange 
are relatively smal 1 (ref. 3). Typical mean free paths for inelastic colli- 
sions are on the order of lO^-lO? cm. 

One plausible explanation for the exponential profiles is that the plasma 
distribution is a direct consequence of the highly viscous flow generated in 
the arcjet nozzle. For comparison, it is worthwhile to examine an alternate 
plume model. The semi-empirical method of G.A. Simons is a widely used tech- 
nique for the prediction of a rocket plume structure (ref. 12). The method 
accounts for the expansion of the viscous, supersonic boundary layer in addi- 
tion to the inviscid core of the nozzle. It is an analogous technique to that 
of Reynolds (ref. 6) except that the major plume parameters are written in 
terms of conditions at the nozzle exit plane and boundary layer thickness 
rather than relying on experimental data obtained in the plume. Simons deter- 
mined that the appropriate selections for the function f(0) were: (1) a 

cosine distribution for the flow in the plume which originates from the invis- 
cid core, and (2) an exponential decay for flow which originates from the vis- 
cous, supersonic boundary layer. No provision was made for subsonic boundary 
layer expansion. The low power arcjet thruster is characterized by highly vis- 
cous nozzle flow. Under typical operating conditions, the Reynolds number, 
based on throat diameter, is <500. One consequence of low Re nozzle flow is 
that the boundary layer makes up a large portion of the nozzle volume and, in 
some cases, may engulf the inviscid core. As the analysis of Simons indicates, 
viscous flow generated in the nozzle decreases less rapidly in the plume than 
does core flow, providing an explanation for the exponential decay in arcjet 
plasma density. 


Plume Impacts Assessment 

Figure 6 displays the relative importance of the three major sources of 
signal attenuation for a 4 GHz signal transmitted through a typical arcjet 
plasma. For the assumed plasma conditions, the refraction loss dominates 
except for densities approaching the critical number density. The reflection 
loss, shown in the figure as a function of n e /n c - only, is the absolute worst 
case, i.e., there is no provision for a slab thickness as in the impacts calcu- 
lations to follow. In practical application where a slab dimension D has been 
defined, the transmission ratio is an oscillatory function which varies between 
(2p/( 1 + h 2 ))2 and 1 as the relative phasing of internal reflections changes. 
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It should be noted that dissipative losses remain low even at high n e /n c 
ratios. The effective collision frequency (assumed here to.be 'the electron-ion 
collision frequency) is on the order of 10^/sec which is an upper bound for a 
1 kW arcjet plume. Because the dissipative, losses are relatively small, absorp- 
tion can be neglected for the situation proposed. 

Figure 7 describes the arcjet plume structure based on the model input 
conditions listed in table I. Plotted in the figure are lines of constant ' 
electron number density emanating from the arcjet source'. The unusual contour 
shape at large angles from centerline is a result of the exponential shape 
function. As mentioned previously, for z < 10 cm, the calculated values of 
electron number density are higher than experimental observations. 

The calculated electron number density that a transmission signal would 
encounter is shown in figure 8 for two propagation paths. For paths which 
intersect the plume centerline, as illustrated by <p Q = 70°, the profile is 
highly collimated with the maximum number density near the center'of the 
assumed plasma slab diameter D. For propagation paths which run parallel to 
the plume centerline (<p 0 = 0°) the maximum number density occurs at the plume 

boundary. The change in the refractive index over the distance of a wavelength 
is the primary cause of reflection at the vacuum-plasma boundary. The sharp 
rise in electron number density at the intersection of the plume boundary is a 
functional limitation of the modelling technique and is probably not real. 
However, it magnifies the reflection problem because, for propagation paths 
which pass near the arcjet source, the higher number density at the boundary 
implies significant reflection of the signal. Table III lists the necessary 
criteria for the validity of the RF signal-plasma interaction model. Sample 
calculations are provided for the above propagation paths. 

- Figures 9 and 10 present the calculated transmission impacts .in terms of 
signal attenuation and phase shift as a function of propagation paths, running 
both parallel to and intersecting the plume axis. The signal loss due to . 
refraction (figs. 9(a) and 10(a)) is presented in units of dB where the trans- 
mission loss is the ratio of power received with the plasma present to the 
power received without the plasma. Except for propagation paths which pass 
very near the arcjet source, the transmission loss is negligible. For propaga- 
tion paths which pass within 15 cm of the source, however, there is significant 
attenuation of the transmitted signal' due' to. refraction of. the beam. Although 
not indicated in the figure, reflection of the signal at the sharp plasma-vacuum 
boundary may also occur for paths close to the source, particularly for paral- 
lel propagation paths. 

The calculated signal attenuation is relatively insensitive to the value 
of slab diameter chosen although, in general , the refraction loss increases 
with decreasing slab diameter. As an example, for the propagation path 
tp Q = 0° , the transmission loss for an assumed slab dimension D. of 542 cm is 

-2x1 0 -4 dB. For 0 = 50 cm, the loss remains negligible at -6xl0~ 3 dB. 

Figures 9(b) and 10(b) present the phase shift as a function of propaga- 
tion path. The shift is a phase advance as the plasma index of refraction is 
less than unity. In contrast to transmi ssion * losses , the phase advance 
increases with increasing signal .frequency. Although there is a noticeable 
phase advance for all propagation paths, the shift is again most severe for 
paths which pass close to-the thruster source. In the antenna far field, a 
phase shift of any magnitude is not significant unless only a portion of the 


11 



plane wave M s shifted in-phase rather than the entire wavefront. Nonuniform 
phase shifts of this sort would limit receiver tracking capability: 

As mentioned previously, the antenna dimensions and the thruster-antenna 
placement were selected such that the RF signal-plasma interaction occurred in 
the antenna far field. In actual flight application this assumption might not 
always be valid as the antenna dimensions may be on the order of 1 m or more. 
In the antenna near field, a complex numerical analysis is required to fully 
define the effects of RF signal-plasma interactions. For example, a phase 
shift of even a few degrees in the antenna near field could effectively dis- 
tort the antenna reflector geometry with severe communications impacts. 


CONCLUDING REMARKS 

This paper has presented a first order evaluation of the potential inter- 
action of a low plasma density arcjet plume with RF signals in common use on 
communications satellites. Modelling the plume as a plasma slab, the communi- 
cations impact for a 4 GHz signal was evaluated for propagation paths running 
both parallel to and intersecting the plume centerline. The impacts were pre- 
sented in terms of signal attenuation and phase shift from transmitting to 
receiving antennae. The dominant transmission loss mechanism was refraction 
of the signal rather than absorption or reflection at the plasma-vacuum bound- 
ary. Unless the propagation path passes very close to the thruster source, 
there appears to be neglibile attenuation due to refraction. A phase advance 
was predicted for all propagation paths. However, in the antenna far field, a 
phase shift of any magnitude is not significant unless only a portion of the 
plane wave is shifted in phase rather than the entire wavefront. 

The two parameter, source flow equation used in .this analysis is adequate 
to conduct an initial assessment of the arcjet plume communications impacts. 
However, a thorough understanding of . the dominant gas/plasma dynamic processes 
occurring in the arc discharge and nozzle regions is required before improve- 
ments to the plume model can be made.. Also, an accurate prediction of plasma 
number density near the plume boundary is needed to fully evaluate the magni- 
tude of signal reflection losses. It is suspected that a more gradual plasma 
boundary exists but measurements in a facility of higher vacuum are needed to 
provide a better description of the complete expansion. Further, this analysis 
investigated only the far field interaction of a, uniform plasma slab with an RF 
signal. A more rigorous calculation of the communication impacts would include 
the effects of a nonuniform plasma density distribution as well as the interac- 
tion processes which occur in the antenna near field. Second order impacts 
such as re-polarization of an EM wavefront would also be considered. 

Lastly, future communications satellites will utilize higher frequency 
bands (20 and 30 GHz) for transmission. This technology advance reduces the 
impacts of’ low power arcjet thruster plumes, because the signal attenuation 
impacts decrease with increasing signal frequency. However, the arcjet plume 
electron number density and temperature increases significantly with power 
level. As the arcjet thruster technology is scaled to increasing power levels 
to meet other space propulsion applications, it cannot be assumed that the two 
technology advances will offset each other in terms of plume' impacts. In fact, 
it seems likely that dissipative effects due to collisions (absorption) will 
also have to be considered in the impacts evaluation of arcjet thrusters above 
a few kW. 
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TABLE I. - ARCJET PLUME MODEL PARAMETERS AND THRUSTER-ANTENNA 
CONFIGURATION INFORMATION FOR IMPACTS EVALUATION 


Arcjet Plume Model Parameters 


Source density, n e> cor-* 3xlo'3 

Constant, A(N,X) 1.1 

Plume parameter, X 5 

Plume parameter, N 0.8 


Antenna Dimensions and Signal Information 


Diameter, D a , cm 40 

Length, L, cm 30 

Signal frequency, u 2.513xl0 10 rads/sec (4.0 GHz) 

Signal wavelength, X s , cm 7.5 

Critical number density, n c , cnr^ 1.984x10" 


Thruster-Antenna Configuration 


Distance behind exit plane, A, cm 50 

Case study 1 - variable <p Q : 

B, cm 200 

<P 0 , deg 0 to 76 

Case study 2 - variable B: 

8, cm 0 to 200 

<f> 0 . deg ^ 0 
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TABLE II. - ASSUMPTIONS NECESSARY FOR APPLICATION OF RF SIGNAL-PLASMA INTERACTION 
MODEL. SAMPLE CALCULATIONS ARE PRESENTED FOR TWO PROPAGATION PATHS. 


[A = 50 cm; B = 200 cm. ] 


Assumption 

Crl teria 

TO = 0° 

To = 

Only free electrons 
interact with the RF 

nig 

CO >> cog. 

uqp = 2.5x10'° rads/sec 
log) = 9x10' 

upp = 2.5xl0 10 rads/sec 
uBi = 9x10' 

signal 

“Pe >> W Pl 

up e = 7x10® 
up} = .5x10® 

upe = 7x10® 
up} = 5x10® 

No absorption loss to 
col 1 i sions 

CO >> v e i 

u R p -= 2 . 5x10'° 
v e } = 2xl0 4 /sec 

vg } = 2xl0®/sec 

URP = 2 . 5x 1 0 10 

Isotropic plasma 

Linear plasma; 1 . e , RF 
power density in the 
plasma low so that 
plasma i s not affected 
by RF wave 

<0 >> «8 e 
Fixed 

assumption 

copp = 2.5xl0'0 
uge = 2x10® 

corp = 2 x 1 0 1 0 rads/sec 
cogg a 2x1 0® 



Adiabatic plasma 

1 Ak 4ir 

k Ax << X s 

r->’ 



. 

i g - 

i g - 



(note: signal loss by 

(signal loss by eq. (6) 

Permeability of plasma = 

D1 amagnetl c 
effect on 
permeability 
Is small 

eq. (6) Is T - 0 dB> 

1sT~ 8x10” 3 dB) 



No induction effects 

D /X >> 1 

d S 

D » 40 cm 

d 

D a » 40 cm 



X $ « 7.5 cm 

X $ ■ 7.5 cm 

Valid use of geometric 
optics limit 

D /X >> 1 

d S 

D = 542 cm 
X s = 7.5 cm 

D = 78 cm 
X $ a 7.5 cm 

Cylinder equivalent to 
a slab of thickness D 

P/D « 1 

D = 542 cm 
p - 6 cm 

D = 78 cm 
p - 14 cm 

Antenna far field 


91 = 375 cm 

gt= 213 cm 



D a = 40 cm 

D a a 40 cm 



X $ = 7.5-cm 

X $ = 7.5 cm 


TABLE III. - ARCJET THRUSTER OPERATING CONDITIONS AND ASSUMED 
PLASMA PARAMETERS FOR COMPARISON OF SOURCE FLOW MODEL 
WITH EXPERIMENTAL DATA 

Arcjet thruster operating conditions 


Arc current, A 10 

Arc vol tage, v ......... 105 

Arc power, W ......... I ...... '. 1050 

Mass flow rate (1:2 mixture of N 2 :H 2 by volume) .■ . . 45 mg/sec 


Assumed plume parameters for model comparlslon 


Source density (at nozzle exit plane) I.lxl0 ,3 /cm® 

A<N,X) 0.87 

Parameter, X 4.0 

Parameter, N 0.8 






FIGURE 1. - ILLUSTRATION OF POTENTIAL ARCJET PLUME INTERACTION 
WITH THE TRANSMISSION OF COMMUNICATIONS SIGNALS. 


r - UNREFRACTED RAY 



FIGURE 2. - GEOMETRICAL OPTICS APPROXIMATION FOR THE EFFECTS OF 
REFRACTION (ADAPTED FROM FIG. R.IO, REF. 7). 
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FIGURE if. - ILLUSTRATION OF THRUSTER-SPACECRAFT CONFIGURATION FOR 
CALCULATIONS OF TRANSMISSION IMPACTS TO MICROWAVE SIGNALS DUE TO 
ARC JET PLASMA PLUME. • 



RATIO n 0 /n 



(B) AXIAL STATION z = 18. A CM. 



(C) AXIAL STATION z = 8.3 CM. 


FIGURE 5. - COMPARISON OF EXPERIMENTAL DATA WITH AN- 
ALYTICAL MODEL. EXPERIMENTAL PROFILES ARE BASED ON 
MEASURED ELECTRON CURRENT NORMALIZED TO CENTERLINE 
VALUE. 
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TRANSMISSION LOSS. dB 


27TUD 7T 



0 

RATIO OF Hq/Hq 


FIGURE 6. - RELATIVE IMPORTANCE OF SOURCES OF TRANS- 
MISSION LOSS AS A FUNCTION OF n e /n c FOR A PLASMA 
TYPICAL OF AN ARCJET PLUME (V = *JX10 6 /SEC; D = 375 CM). 


DISTANCE ALONG PLUME CENTERLINE/ CM 

FIGURE 7. - CONSTANT NUMBER DENSITY CONTOURS. SOURCE 
DENSITY = 3x10 13 /CM 3 ; N = 0.8; X = 5. 





PHASE SHIFT, DEG TRANSMISSION LOSS, dB 



TION ANGLE <P Q . AN ANGLE OF <p Q = + 75° CORRESPONDS 
TO A PROPAGATION PATH DIRECTLY INTO THE ARCJET 
SOURCE. 



(B) PHASE SHIFT AS FUNCTION OF PROPAGATION ANGLE <p Q . 
FIGURE 9. - TRANSMISSION IMPACTS FOR PROPAGATION PATHS 
INTERSECTING PLUME CENTERLINE. [A = 50 CM; B = 200 CM]. 
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0 40 80 120 160 200 

DISTANCE B, CM 

(B) PHASE SHIFT AS A, FUNCTION OF DISTANCE B. 

FIGURE 10. - TRANSMISSION IMPACTS FOR PROPAGATION PATHS 
PARALLEL TO PLUME CENTERLINE. [A - 50 CM; (p Q = 0°]. 
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